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ABSTRACT 

We model populations of isolated neutron stars in the Galaxy following their orbital and 
magneto-rotational evolution. It is shown that accretors become more abundant than coolers 
at fluxes below '^ 10~^^ erg cm~^ s~^, and one can predict that about one accreting neutron 
star per square degree should be observed at the Chandra and Newton flux limits of ~ 10~^^ 
erg cm~^ s~^. The soft ROSAT sources associated with isolated neutron stars can be relatively 
young cooling objects only if the neutron star birth rate in the Solar vicinity during the last 
~ 10^ yr is higher than that inferred from radiopulsar observations. 

Subject headings: accretion, accretion disks — stars: kinematics — stars: magnetic fields — 
stars: neutron — stars: statistics — X-rays: stars 



1. INTRODUCTION 

Despite intensive observational campaigns, no ir- 
refutable identification of an isolated accreting neu- 
tron star (NS) has been presented so far. Six soft 
sources have been found in ROSAT fields which are 
most probably associated to isolated radioquiet NSs. 
Present X-ray and optical data however do not allow 
an unambiguous identification of the physical mech- 
anism responsible for their emission. These sources 
can be powered either by accretion of the interstellar 
gas onto old (~ 10^'' yr) NSs or by the release of inter- 
nal energy in relatively young (~ 10^ yr) cooling NSs 
(see Treves et al. 2000 and Motch 2000 for recent re- 
views). The ROSAT candidates, although relatively 
bright (up to ~ 1 cts~^), are intrinsically dim and 
their inferred luminosity (L ~ 10'^^ ergs~^) is near 
to that expected from either a close-by cooling NS 
or from an accreting NS among the most luminous. 
Their X-ray spectrum is soft and thermal, again as 
predicted for both accretors and coolers (Zane, Tur- 
olla & Treves 2000; Treves et al. 2000). Up to now 
only two optical counterparts have been identified 
(RXJ 1856, Walter & Matthews 1997, for which a 
distance estimate of ^ 60 pc has been very recently 
obtained, Walter 2000, and RXJ 0720, Kulkarni & 
Van Kerkwick 1998). In both cases an optical excess 
over the low-frequency tail of the black body X-ray 
spectrum has been reported. While detailed multi- 
wavelength observations with next-generation instru- 



ments may indeed be the key for assessing the true 
nature of these sources, other, indirect, approaches 
may be used to discriminate in favor of one of the 
two scenarios proposed so far. 

In this Letter we compute and compare the log A^ - 
log S distribution of both accreting and cooling NSs, 
to establish the relative contribution of the two pop- 
ulations to the observed number counts. Previous 
studies derived the log N - log S distribution of ac- 
cretors (Treves & Colpi 1991; Madau and Blaes 1994; 
Manning et al. 1996) assuming a NSs velocity distri- 
bution rich in slow stars {v <^ 100 kms^ ). More 
recent measurements of pulsar velocities (e.g. Lyne 
&: Lorimer 1994; Hansen & Phinney 1997) and upper 
limits on the observed number of accretors in ROSAT 
surveys (Danner 1998a, b) point, however, to a larger 
NS mean velocity (see Treves et al. 2000 for a critical 
discussion). Recently Neiihauser & Triimper (1999, 
NT99 hereafter) compared the number count distri- 
bution of the ROSAT isolated NS candidates with 
those of accretors and coolers. Here we address these 
issues in greater detail, also in the light of the lat- 
est contributions to the modeling of the evolution of 
Galactic NSs (Popov et al. 2000, P2000 hereafter). 
In §2 we use a population synthesis model, devel- 
oped in P2000, to trace the properties of the NS spin, 
magnetic field and accretion luminosity in the Galac- 
tic potential using a detailed map of the interstellar 
medium (ISM). Cooling NSs are explored separately 
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in §3 within a simpler model of local sources. Results 
are discussed in §4. 

2. ACCRETING ISOLATED NEUTRON STARS 

The census model adopted here follows closely that 
developed in P2000. All NSs are born in the Galac- 
tic plane, their birth rate is taken to be proportional 
to the square of the ISM density and is constant in 
time. Initially a NS has a circular velocity fixed by 
its position in the Galactic plane and determined by 
the Galactic potential. To this, an additional kick ve- 
locity with random orientation and modulus selected 
from a Maxwellian distribution is imparted. 

Time resolution is adapted to follow in detail the 
different evolutionary phases the star experiences. 
NSs are born in the ejector (E) phase and then, as 
braking due to electromagnetic torques increases the 
period, come to the propeller (P) phase. Further 
spin-down, produced by accretion torques, may drive 
the star in the accretion (A) phase when matter can 
penetrate the magnetosphere reaching the stellar sur- 
face (see e.g. Lipunov 1992). Accretion is then as- 
sumed to proceed at the local Bondi rate and the 
large scale ISM distribution is mapped as in P2000. 
The small scale structure of the ISM (w lOOAU-1 pc, 
see e.g. Meyer k, Lauroesh 1999; Treves et al. 2000), 
has not been taken into account as it would produce 
an unacceptable increase of the computing time, lead- 
ing to a worse statistics. A detailed mapping of the 
ISM on all scales is fundamental in assessing the ob- 
servability of single, luminous sources (see Zane et 
al. 1996), but the effects produced by its multiphase 
structure on the evolution of the entire NS population 
are likely to average out. 

Here we consider only the simplest picture, where 
the star magnetic field is constant and the dipole mo- 
ment ^ has a log-gaussian distribution, similar to that 
observed in radiopusars ((log/i) = 30.06, a = 0.32 in 
CGS units). 

The magneto-rotational evolution is followed using 
essentially the same approach as in P2000. In partic- 
ular, the typical values of the transition periods are 
~ 25 s and ~ 1500 s for E^ P and P^ A transitions, 
respectively. As soon as the star enters the accretor 
stage, the period is set to the value typical of an ac- 
cretor embedded in a turbulent medium (Lipunov & 
Popov 1995; Konenkov k Popov 1997). 

The statistics of accretors was improved by a factor 
> 10 relative to P2000 by evolving only stars from the 
low- velocity tail {v < 100 km s~^) of a Maxwellian 
with (f) = 300 km s~^. This assumption appears 
justified since the time spent in the ejector phase in- 
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ing NSs (P2000). Results can be easily scaled to any 
{v) > 200 km s-^ 

Nearly all highly magnetized isolated NSs in the 
velocity range considered here become accretors, and 
in calculating the log A'^ - log S distribution we ob- 
tain good statistics for fluxes < 10^^^ erg cm^^ s^^. 
At larger fluxes the number of sources is too small 
to give reliable counts, but the extrapolation log A^ ^^ 
— (3/2) log S can be safely used for log S" ;^ —11. 

Our main results are presented in Figures 1 and 3, 
and refer to a total of 10^ NSs in the Galaxy. The 
logA^ - logS" distribution of accretors is computed 
for sources within a distance of 5 kpc, together with 
the distributions of velocity, accretion rate and ef- 
fective temperature (calculated adopting as effective 
surface the polar cap area). Isolated NSs with fields 
< 0.5 X 10^^ G never appear as accretors in any run. 
The brightest accretors have luminosities ^ 10^^ erg 
s~^, but for the majority of them L ~ 10^^ — lO^'^ erg 
s~^. Our sources are relatively hot, with a mean ef- 
fective temperature 300-400 eV and for them absorp- 
tion can be neglected. We predict, on average, about 
1 source per square degree, in the energy range 0.5-2 
keV, for limiting fluxes about 10~^^-10~^^ erg cm~^ 
s^^. This implies --^ 3 x 10^ sources in the whole sky. 
Note that they are significantly concentrated toward 
the Galactic plane with a center-anticenter asymme- 
try. 
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Fig. 1. — Upper left panel: the logiV - log 5 distribu- 
tion for accretors within 5 kpc from the Sun. The two 
curves refer to emission from the entire star surface and 
to polar cap emission in the range 0.5-2 keV; two straight 
lines with slopes -1 and -3/2 are also shown for compari- 
son. From top right to bottom right: the velocity, effective 
temperature and accretion rate distributions of accretors; 
all distributions are normalized to their maximum value. 



present choice for (v) follows from radiopulsar obser- 
vations and current limits on the detection of accret- 
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3. COOLING NEUTRON STARS 

Soon after their birth in a type II supernova event, 
cooHng NSs emit X-ray luminosities as high as ~ 10^^ 
erg s~^ for a few thousands years. Before disappear- 
ing from the X-ray sky, they shine as faint soft sources 
(L '-^ 10^^ erg s~^) for about a mihion year (e.g. 
Yakovlev et al. 1999). 

In the determination of the logA^-logS*, the cooler 
sample should in principle be drawn from the same 
parent population as accretors. However, as their to- 
tal number is quite low (~ 10~^ of the total number 
of NS in the Galaxy) owing to the short duration of 
the cooling phase, the results would be affected by 
the very poor statistics. Notice however that coolers 
can travel only a distance ~ 400 pc during their life- 
time and have softer spectra than accretors, so they 
appear mostly as a local population of sources (due 
to interstellar absorption). 

For this reason, we prefer to approximate their spa- 
tial distribution around the sun as homogeneous with 
a scale height of ^^ 450 pc. We take cooling NSs as 
"standard candles" with L = 10^^ erg s^^, giving 
off a black-body spectrum with effective temperature 
^ 50 eV. The duration of the cooling phase is taken 
to be 10^ yr, according to the "slow cooling" sce- 
nario (see NT99). We used two typical values for 
the total NS spatial density: 0.33 x 10^^ pc~^ and 
3.3 X 10~^ pc~^. The former estimate follows from ra- 
diopulsar statistics and was used by NT99, while the 
latter corresponds to a total population of 10^ NSs, as 
implied by constraints on nucleosynthetic yields and 
historic supernova rates (Arnett et al. 1989; van den 
Bergh &: Tammann 1991). Accordingly, coolers have 
a mean density of 0.33 — 3.3 x 10~^pc~^. 

The ISM structure is treated in a very simple way. 
The ISM density inside the Local Bubble of radius 
ri = 140 pc (Sfeir et al. 1999) is n = 0.07 cm'^. At 
larger distances the interstellar gas is assumed to be 
uniformly distributed in a disk of half-thickness 450 
pc with n = 1 cm~^. 

Results are shown in Figure 2, where we compare 
the observed distribution of ROSAT candidates with 
that of coolers, as predicted by the present model (in- 
cluding absorption) and by NT99, the latter obtained 
for a total NS number Ntot = 2 x 10^. As it is appar- 
ent, coolers dominate over accretors at large fluxes 
(for the same Ntot), essentially because of their much 
higher luminosity (10^^ vs. 10^^ — lO^'^ erg s~^). 

This simple model reproduces the most important 
feature, the "flattening" of the log N - log S distribu- 
tion outside the Local Bubble. Such strong flattening 
can help to explain the observed data, if one assumes 
that most of the brightest isolated NSs are identified 
(Schwope et al. 1999). 

The value of the flux at which the distribution flat- 



tens out (the "knee") depends on the size of the Lo- 
cal Bubble r/. The ISM densities and the duration of 
the cooling time affect instead the number of sources. 
Since the cooling time can be an important parame- 
ter, the comparison of theoretical and observed dis- 
tributions could, in principle, help in discriminating 
between "fast" and "slow" cooling models. 
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Fig. 2. — The logA^-logS' for coolers. The upper curve 
(filled circles) correspond to a NS density of 3.3 x 10^'^ 
pc^"^ and the lower one (open circles) to 3.3 x 10~^ pc 
here ri = 140 pc. 
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4. DISCUSSION 

The statistical analysis presented in the previous 
sections shows that both scenarios, accretors and 
coolers, have difficulties in explaining, under usual 
assumptions, the observed properties of isolated NS 
ROSAT candidates. The main problem is to pro- 
duce enough sources that are, at the same time: i) 
relatively bright (high count rates, ^ 0.1 cts s~^); 
ii) close (low absorption Njj ^ 10^'^cm^^); iii) soft 
{Teff ~ 50 — 100 eV); iv) slowly rotating (for two 
candidates, RX J0720 and RX J0420, the detected 
periods are 8.4 and 22.7 s). 

Polar cap accretion for a non-decaying magnetic 
field, B Ri 10^^ G, cannot produce enough bright 
sources and the emitted spectrum is likely to be 
harder, see Figure 1. Moreover, accretors are ex- 
pected to have much larger periods, P ^ 10^ s. 
On the other hand, relatively young cooling NSs 
can be spun down to periods ~ 20 s in < 10^ yrs 
(when they are still hot) only for B ^ 10^^ G, since 

P ~ 15 Bu{t/ 10'^ yr)i/2 s for dipole losses. If ah NSs 
experience an initial radiopulsar phase, their total 
number in the Galaxy would be about (1 — 3) x 10^ 
(Lyne et al. 1998). In this case there will be not 
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bright enough isolated NSs in the Solar vicinity in 
both scenarios. 

If Ntot ~ 10^, coolers can explain the large number 
of bright close-by sources and their typical temper- 
atures (see Figure 2). However, the problem of re- 
producing the observed periods still remains, unless 
the progenitors of the ROSAT pulsating candidates 
are NSs born with long periods, which never passed 
through the pulsar phase; alternatively they could 
be magnetars with ultra-high fields, as suggested by 
Heyl & Kulkarni (1998). If proved correct, the "fast" 
cooling scenario, in which the duration of the active 
X-ray phase is shorter, would make again problematic 
to account for the observed bright sources. 

In the accretion scenario, increasing the total num- 
ber of NSs is also a need. As already shown by Livio 
et al. (1998) and Colpi et al. (1998), magnetic field 
decay can explain the low temperatures and periods. 
Decay can both increase and decrease number of ac- 
creting isolated NSs (Popov, & Prokhorov 2000). But 
a complete statistical analysis is still to come. 
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Fig. 3. — Comparison of the log N - log S distri- 
butions for accretors and coolers together with observa- 
tional points, the naive logA^ - logS* from P2000 and the 
ROSAT Bright Survey (RBS) limit (Schwope et al. 1999). 
The scale on the top horizontal axes gives the flux in erg 



In any case, using "standard" assumptions on the 
velocity, spin period and magnetic field parameters, 
the accretion scenario can not explain the observed 
properties of the six ROSAT candidates. 

Note that the number of accretors (of all luminosi- 
ties) in the Galaxy is found to be about two orders 
of magnitude larger than the number of young cool- 
ers: few percents versus 0.01% of the total number 
of isolated NSs. But on average accretors are 3 or- 
ders of magnitude fainter than a typical cooler. A 
key result of our statistical analysis is that accretors 
should eventually become more abundant than cool- 
ers at fluxes below 10~^^ erg cm~^ s~^. This is illus- 
trated in Figure 3 in which we summarize the main 
outcome of our simulations. 

ACKNOWLEDGMENTS 

We wish to thank J. Lattimer, E. van den Heuvel, 
V.M. Lipunov, and S. Campana for useful discus- 
sions. SBP and MEP also thank the University of 
Insubria for financial support and the Universities 
of Milano-Bicocca, Padova and the Brera Observa- 
tory (Merate) for their kind hospitality. The work of 
SBP and MEP was supported through grant RFBR 
00-02-17164 and NTP Astronomy grants 1.4.4.1. 
and 1.4.2.3. Support from the European Commis- 
sion under contract ERBFMRXCT98-0195 and from 
MURST under contract COFIN98021541 is acknowl- 
edged. 



10' 



REFERENCES 



Arnett, W.D., Schramm, D.N. & Truran, J.W. 1989, ApJ, 339, 

L25 
Colpi, M., TuroUa, R., Zane, S., & Treves, A. 1998, ApJ, 501, 

252 
Danner, R. 1998a, A&AS, 128, 331 
Danner, R. 1998b, A&AS, 128, 349 

Hansen, B.M.S., & Phinney, E.S. 1997, MNRAS, 291,569 
Heyl, J.S., & Kulkarni, S.R. 1998, ApJ, 506, L61 
Konenkov, D.Yu., & Popov, S.B. 1997, AZh, 23, 569 
Kulkarni, S.R., & van Kerkwijk, M.H. 1998, ApJ, 507, L49 



Lipunov, V.M. 1992, Astrophysics of Neutron Stars (Berlin: 

Springer & Verlag) 
Lipunov, V.M., & Popov, S.B. 1995, AZh, 71, 711 
Livio, M., Xu, C, & Frank, J. 1998 ApJ, 492, 298 
Lyne, A.G., & Lorimer, D.R. 1994, Nature, 369, 127 
Lyne, A.G., et al. 1998, MNRAS, 295, 743 
M:adau, P., & Blaes, O. 1994, ApJ, 423, 748 
Manning, R.A., Jeffries, R.D., & Willmore A.P. 1996, MNRAS, 

278, 577 
Meyer, D.M., & Lauroesh, T.J. 1999, ApJ, 520, L103 



Popov et al. 



Motch, C. 2000, to be published in the Proceedings of the 
Meeting "X-ray Astronomy 1999: Stellar Endpoints, AGN 
and the DiflFuse Background" (astro-ph/0008485) 
Neiihauser, R., & Triimper, J.E. 1999, A&A, 343, 151 (NT99) 
Popov, S.B., & Prokhorov, M.E. 2000, A&A, 357, 164 
Popov, S.B., Colpi, M., Treves, A., TuroUa, R., Lipunov, V.M., 

& Prokhorov, M.E. 2000, ApJ, 530, 896 (P2000) 
Schwope, A.D., Hasinger, G., Schwarz, R., Haborl, F., & 

Schmidt, M. 1999, A&A, 341, L51 
Sfeir, D.M., Lallemont, R., Crifo, F., & Welsh, B.Y. 1999, A&A, 
346, 785 



Treves, A., & Colpi, M. 1991, A&A, 241, 107 

Treves, A., TuroUa, R., Zane, S., & Colpi, M. 2000, PASP, 112, 

297 
van den Bergh, S., & Tammann, G.A. 1991, ARA&A, 29, 363 
Walter, F., & Matthews, L.D. 1997, Nature, 389, 358 
Walter, F. 2000, ApJ, in the press (astro-ph/0009031) 
Yakovlev, D.G., Levenfish, K.P., & Shibanov, Yu.A. 1999, Phys. 

Usp. 42, 737 
Zane, S., Zampieri, L., TuroUa, R., & Treves, A. 1996, A&A, 

309, 469 
Zane, S., TuroUa, R., & Treves, A. 2000, ApJ, 537, 387 



